Simulating quantum operations with mixed environments 
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We study the physical resources required to implement general quantum operations, and provide 
new bounds on the minimum possible size which an environment must be in order to perform certain 
quantum operations. We prove that contrary to a previous conjecture, not all quantum operations 
on a single-qubit can be implemented with a single-qubit environment, even if that environment is 
initially prepared in a mixed state. We show that a mixed single-qutrit environment is sufficient to 
implement of a special class of operations, the generalized depolarizing channels. 
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Future quantum computers may be useful in studying 
the behavior of open quantum systems and the nature 
of decoherence [^^. Instead of performing real experi- 
ments on quantum systems, a single quantum computer 
can be used as an efficient, multiple-purpose simulator for 
a wide variety of physical systems. In general, an impor- 
tant goal of such investigations will be to understand the 
effects arising from interactions between the system of 
interest S and another quantum system E. For example, 
a quantum computer can be used to simulate quantum 
systems S in thermal equilibrium |^ , but such a simula- 
tion requires an additional quantum system E, coupled 
in a particular way to S, to mimic the thermal bath of 
the system. In other applications involving the simu- 
lation of nonequilibrium quantum properties, S could, 
for a molecules whose isomerization dynamics we wish 
to study, represent the relevant conformational states, 
which couple to other molecules E through long-range 
electronic dipolar interactions. In all these applications, 
we wish to to implement E with the smallest quantum 
resources possible, and this Letter investigates the most 
efficient implementation of such quantum environments. 

Suppose S exists in a Hilbert space 7i„ of dimension 
n, and E is in Tir of dimension r. It is well known that 
any quantum operation Q on 7Y„, resulting from some 
interaction with E in 7i,. with arbitrary r, can be per- 
formed by appending a state in 7i„2 , evolving unitarily, 
and then tracing over 7i„2 . The difference between r and 
can represent a significant reduction, since E can be 
a large bath (for example, of harmonic oscillators), and 
r can be infinite. 

Can a general quantum operation be implemented with 
an environment even smaller than dimensions? Lloyd 
conjectured |^ that it is possible to implement a general 
quantum operation on k quantum bits (qubits) with a 
/c-qubit environment - if one prepares the environment 
not in a pure state, as is usually the case, but rather in 
an arbitrary mixed state. 

Here we provide a specific counterexample to this con- 
jecture for fc = 1, although we find that at least for some 



operations, fewer resources are required than was pre- 
viously known. Our counterexample is part of a class 
known as the generalized depolarizing channels, for which 
we show that a three-dimensional environment is suffi- 
cient for simulation. The proof of the counterexample is 
established by the technique of computing Grobner bases. 

Our results also address the following question: sup- 
pose a physical system is given as a black box ~ we can 
prepare S in an arbitrary initial state, and then mea- 
sure the final state of S after a fixed evolution period. 
What is the largest environment E with which S might 
have interacted in this system? A method to completely 
determine the quantum operation x performed by this 
system is known This work goes one step further, by 
showing a way to turn knowledge about x ij^to bounds 
on the nature of E. 

We begin by summarizing the mathematical formalism 
of quantum operations. The most general transformation 
allowed by quantum mechanics for an initially isolated 
quantum system is a linear, trace-preserving, completely 
positive map. Such a map x '-^n — > -A.m, where An 
is the set of operators on a Hilbert space Hn, can be 
decomposed into a set of at most nm m x n matrices Ai 
H (which we shall refer to as "operation elements" ) as 



(1) 
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n and m are the dimensions of the input and output 
Hilbert spaces, respectively. The trace-preserving prop- 
erty implies that the Ai obey the constraint 



(2) 



with In the identity matrix on 7i„. Following Choi 
1^ , the set of all such maps x ■ — > Am we call 
TCP[n, m]. A physical implementation of these maps 
is represented in Fig. ^: A unitary operation on the state 
p (8) |0)(0| (where |0) represents some pure state in an 
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m^-dimensional environment) is performed and then nm 
"degrees of freedom" are traced out: 



x(p) = ^(e..|C/ 



ek) 



(3) 



fc=i 



Here {|efe)}^™i is a set of basis vectors for Tinm- As there 
are at most nm operation elements, it follows that one 
can implement any map in TCP[n, m] with an environ- 
ment of dimension m^. To determine the dimension of 
the parameter space of TCP[n, m] we note that the map 
X does not uniquely determine the set {Ai}™^. Any set 
of TO X n matrices {i3i}™| and {Aj}™ that are related 
by a unitary transformation 



nm 



(4) 



implement the same map x- This freedom corresponds 
to a unitary rotation [/' (see Fig. ^ of the environment 
qubits after the completion of the interaction U. It is 
shown in Q that this unitary equivalence is the only 
freedom in the choice for the set of operators {^ij^ 
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FIG. 1. Implementation of the map x using a pure state 
environment. 

The dimension of the parameter space of all maps 
in TCP [n, to] that can be implemented with a d- 
dimensional environment will therefore be 



7-jn — >m 

purCjrf 



parameters in {A^} 

2^ 



unitary freedom 

[nd/m)^ 



constraint (2) 

(5) 

Thus we have 



where d is such that to divides nd. 

In a more general physical implementation, however, 
the initial state of the environment can be an arbitrary 
density matrix. Consider the set of completely posi- 
tive trace-preserving linear maps x ■ -^n — >■ Am that 
are implemented by an environment that is initially in 
some d-dimensional density matrix. We call this set 
[d, n, in]. The action on the input state p is 



dn/r, 



(6) 
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where {Aj, are now the eigenvalues and eigenvec- 

tors of the mixed environment state. We identify a set 



of TO X n matrices {j4jfc}^!f"{,™j^ in the representation of 
Eq. d): 

Ajk = ^MU\]). (7) 
Unitarity implies that these matrices are constrained, 

Y^^lkAjk^S^Kln- (8) 



There is a residual unitary freedom in choosing the set 

of matrices {Ajk}f^^i"li- The set {Bjm}f=tlm=i '^i*^ 
Bjm = J^k^mk^jki whcrc the c?ri/TO-dimensional uni- 
tary matrix U' does not depend on the label j, imple- 
ments the same quantum operation and also obeys con- 
straint (j^). As before, this freedom corresponds to a uni- 
tary transformation on the environment after the com- 
pletion of the operation. The dimension of the parameter 
space of S'inix[c^, m] can be bounded: 



pure, a — mix.d — purc,d^ 



(9) 



The upper bound is given by the fact that one can always 
simulate a d-dimensional mixed environment with a d^- 
dimensional pure environment. 

From Eq.d) and Eq.® it follows that an environment 
of dimension d < m cannot be used to implement all 
maps in TCP [n, to]. In fact a large set of maps, the 
extremal maps in TCP [n, to], cannot be simulated with 
d < m. A map x that is decomposable in to or fewer 
linearly-independent operation elements is extremal |^ 
in TCP [n, to] . These maps can be implemented with 
a pure-state environment of dimension to; moreover, we 
prove that there does not exist a more efficient implemen- 
tation of these maps using a mixed-state environment: 

Extremality implies that the map x cannot be written 
as a convex combination of linearly independent maps 
X* that each have operation elements {^J} for which 

J2j A'j^ A'j — In for each i. This ensures that only one of 
the eigenvalues in constraint (^) is non-zero, but this in 
fact corresponds to a pure-state environment of dimen- 
sion TO. An example of such an extremal map is a von 
Neumann measurement on a n-dimensional system. The 
set of projection operators {PiY^^i can be implemented 
minimally by using an n-dimensional pure state. 

We now turn to the question of whether all maps in 
TCP[n, to] can be implemented with d = m. Note that 
our parameter count does not exclude this. In the fol- 
lowing, we restrict ourselves to the case n = to = 2. We 
study which maps can be implemented using a single- 
qubit environment and provide a proof that a particular 
qubit channel, the two-Pauli channel, cannot be imple- 
mented in this way. 

We consider a special set of maps, the generalized de- 
polarizing channels 0, which are described by the set 
where 
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(10) 



such that ei + £2 + £3 + £4 = 1 ^i^d the operators Ai 
are given by Ai = I2, A2 = ct^, A^ = (Ty, A4 = ci^. 
One can represent this family of maps geometricaUy as a 
tetrahedron, which is embedded in a cube with vertices 
at (1,-1,-1), (-1,1,-1), (1,1,1) and (-1,-1,1). The 
transformation that relates the parameters ei, £2, £3, £4 to 
the (x, y, z) coordinates is given by x = £1 + €2 — £3 — £4, 
y — £1 — £2 + £3 — £4, and z = £1 — £2 — £3 + £4- The vertices 
of the tetrahedron correspond to a single-operator map. 
Its edges are two-operator maps, the four faces represent 
all three-operator maps, and the points in the interior 
of the tetrahedron are all the four-operator maps of Eq. 

A computer search suggests that only a subset of these 
maps can be simulated by using a qubit environment. 
For this subset we are able to construct an explicit qubit 
solution. At web address |^ one can find pictures of 
the three-dimensional volume that is described by the 
solution set and a picture of the solution set as gener- 
ated by the computer search. The computer work also 
suggests that the dimension of S'niix[2, 2, 2] is equal the 
upper bound of Eq.(^), namely D'^'^^^^ — 12. We find 
this by randomly sampling in the space of all superop- 
erators, that is, we choose random orthonormal vectors 
that make up the columns of the unitary matrix U of 
Eq. (|3|); a finite percentage could be implemented with 
a qubit environment. Thus there is enough "room" for a 
solution, but it is not in the right place, as we will see. 

This solution is constructed in the following way. We 
start with the center of mass of the tetrahedron, the 
point (£1, £2, £3, £4) = (1/4,1/4,1/4,1/4). This channel 
has the property that it maps every input state p onto 
^l2. It can thus be easily implemented by performing 
a SWAP gate on a environment qubit that is initially 
in the ^12 state and the input qubit. The SWAP gate 
on two registers \a)\h) gives |fo)|a). Then one consid- 
ers the line that departs from a vertex, say the point 
(eij £2, £3, £4) — (1,0,0,0), and goes through the center 
of mass. This one-dimensional set of channels is char- 
acterized by £2 = £3 = £4 and represents the regular 
depolarizing channel 0. Performing a V SWAP on a 
il2 environment and the input qubit implements these 
channels, up to £1 = 1/4. The integer m is related to 
the £ parameters by £2 = £3 = £4 = sin^(2^)/4. One 
extra step of generalization gives us an even larger set of 
channels. The unitary matrix is a somewhat generalized 
form of VSWAP, 



U = 



( e''^ cos (/>i ie'^ sin \ 

cos 4>2 i sin 02 

i sin (j)2 cos 02 

\ ie'^ sin 0i e*^ cos 0i / 



, (11) 



and the environment is again prepared in state 512- We 



can determine the operation elements and express these 
as linear (unitary) combinations of the Pauli matrices. 
This leads to an expression of the parameters £i in terms 
of (6*, 01 , 02) e [0, 27r] X [0, 27r] x [0, 2tt] : 



£1 |(C0S2 01 

£2 = \{s\t? 01 

£3 = Ksin^ 01 

£4 = l(cos2 01 



cos 

2 



- 2 cos 01 COS 02 cos ( 
sin~ 02 + 2 sin 0i sin 02 cos 9 
■ sin^ 02 — 2 sin 0i sin 02 cos d 

- COS^ 02 ~ 2 COS 01 cos 02 COS i 



(12) 



We now turn to another set of maps, the two-Pauli 
channel, which is given by the three operators 



Ai 



l2V^,^2 = a^^(\-x)l2,A-i = iay^J{\~x)l2. 

(13) 



We will prove that for < x < 1, there is no qubit 
environment which simulates this channel. For x — or 
X — 1 there is a two-dimensional environment that can 
simulate the channel as the channel has two operation 
elements when x = and only one operator when a: = 1. 

Any unitary linear combination of the Ai , A2 and A3 
may be written as 



Bk = 



hh\fx 
ak 



(cfc - afe).y/i(l - x) 



(14) 



with appropriate constraints resulting from unitarity on 
the coefficients ak,bk,Ck- This new set of operators 
{Bk}\^i will implement the same channel due to Eq.(H). 
Furthermore, these operators Bk are constrained through 
Eq. (^). For notational convenience, we define 
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V2 
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ci - ai 
bo 
bi 



\Wl) = -TS 



1 
V2 



^1 



0-2 

0-3 

C2 - 

C3 - 
62 
63 



C2 
C3 
0-2 
03 



(15) 



Using the assumption 7^ a; 7^ 1 and by linearly com- 
bining all the equations we obtain: 



{vo\wo] 


+ {uo\vo) = 


91 





(16) 


{vi\wi) 


-1- {ui\vi) = 


92 = 





(17) 


{vo\wi) 


-1- {uo\vi) = 


53 = 





(18) 


{uJo\vi) 


-1- {vo\ui) = 


.94 = 





(19) 


{uo\uo) 


- {wo\wq) = 


95 = 





(20) 


(mi|mi) 


- {wi\wi) = 


56 = 





(21) 


"o|uo) + ( 


ui\ui) - 1 = 


57 = 





(22) 


{vo\vq) + 


vi\vi) - 1 = 


58 = 





(23) 


{y-o\vo) 


-1- {ui\vi) = 


59 = 





(24) 


{uo\wo) 


-1- {ui\wi) =: 


510 = 





(25) 


{uo\ui) 


- {wo\wi) =: 


511 = 





(26) 
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Writing each of the coefRcients a^, bk, and Ck in the 
form Xj+ixj+i (where = —1), we get a system of poly- 
nomial equations Re((7i) = Im((7i) = . . . = Im(5ii) — 0, 
where Re((7fc) and Im((7fc) are polynomials in the vari- 
ables a;i,...,X24 with real coefficients. To show that 
this system of equations has no solution we make use 
of Grobner bases (see e.g. Q). The computation of a 
Grobner basis with Buchberger's algorithm generalizes 
the Euclidean algorithm to compute the greatest com- 
mon divisor (GCD) of univariate polynomials pi{x) and 
P2{x). In that case, the GCD g{x) can be written as 
a "linear" combination g(x) = fi(x)pi{x) + f2{x)p2{x). 
The two univariate polynomials pi and p2 have a com- 
mon root if and only if their GCD is non-trivial, i.e., 
gix)^l. 

For multivariate polynomials, a common solution ex- 
ists iff the Grobner basis of the ideal generated by them 
is non-trivial, i.e., does not contain a constant. In our 
case, using the computer algebra system Magma 
we have shown that there exist polynomials /i, . . . , /n 
such that /j (2^1, •■• ,2:24)ffi(a;i, ... ,X24) = 1, i.e., 

the Grobner basis contains 1 and there is no solution of 
the equations (|6l)-(|26|). □ 



Despite the above proof, it turns out that the class 
of channels we have been studying do not require a 
two qubit environment {d = 4) for their simulation; a 
mixed qutrit (d = 3) suffices. For generalized depolar- 
izing channels, there will be nine operators, 
We set one eigenvalue A3 = and thus A31 — A32 — 
^33 = 0. If £162 > £364 the solution is An — 0, 

A21 = \/i-2- ^i^i/^icrx, A12 =j/e^z, A22 = \JTl<Jy^ 

^13 — \/eil2, and A23 = —i\J ^z^aI ^\Ox- Otherwise, we 
take All = 0, A21 = V^4 - eie2/e3crj^, A\2 = \/ei'i-2, 

, and A23 = i^/eie^/e^ay. 



^22 = \f^2^x, ^13 — 

One can check that this set implements any generalized 
depolarizing channel and satisfies Eq. (||). 

On the basis of the computer work we conjecture that 
any map in TCP [2, 2] can be simulated with a qutrit 
environment. Also, the numerics suggest that one can 
always set one eigenvalue to zero. Furthermore, we have 
some numerical evidence that channels that have three 
linearly independent operation elements can never be 
simulated with a qubit environment. 

Our results provide new bounds on the size of an envi- 
ronment needed to simulate certain quantum operations 
on single qubits. However, we have only addressed simple 
mappings on the smallest input space. Many questions 
now arise: how do these results generalize to mappings 
on n-dimensional systems? A relevant scenario might be 
n uses of the generalized depolarizing channel, where the 
environment can be shared between the channels. In such 
a case, might a qubit environment per channel suffice for 
large n? A nice extension of the generalized depolariz- 
ing channels are the channels that are defined with the 
Heisenberg group elements These channels on n- 



dimensional inputs are mixtures of a set of r? unitary 
matrices U(i,j). However, it is not straightforward to 
construct solutions, as in the qutrit case, for a general 
"Heisenberg channel," and we have no insight at the mo- 
ment of what gain one can get by using mixed states here. 
The questions we have formulated also apply to the con- 
struction of generalized measurements: how large an en- 
vironment is needed for the minimal-size construction of 
arbitrary generalized measurements on an n-dimensional 
system? We hope our results and the questions they mo- 
tivate will be useful in future quantum computing appli- 
cations, and provide fundamental insights into the often 
strange properties of quantum systems. 
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